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Bouncing Anisotropic Universes with Varying Constants 
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We examine the evolution of a closed, homogeneous and anisotropic cosmology subject to 
a variation of the fine structure 'constant', a within the context of the theory introduced 
by Bekenstein, Sandvik, Barrow and Magueijo, which generalises Maxwell's equations and 
general relativity. The variation of a permits an effective ghost scalar field, whose negative 
energy density becomes dominant at small length scales, leading to a bouncing cosmology. 
A thermodynamically motivated coupling which describes energy exchange between the ef- 
fective ghost field and the radiation field leads to an expanding, isotropizing sequence of 
bounces. In the absence of entropy production we also find solutions with stable anisotropic 
oscillations around a static universe. 
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I. INTRODUCTION 

Spatially homogeneous cosmological models 
are a key area of study within relativity. The in- 
troduction of anisotropies gives rise to models in 
which a richer dynamical structure emerges, yet 
the cosmology remains simple enough to provide 
analytic and simple numerical results. These 
models serve as a test-bed for physical theories, 
and allow us analyse questions about why the 
universe appears to be highly isotropic, whether 
inflation occurs for generic or stable sets of initial 
data, the effects of anisotropy on astronomical 
observables, and the behaviour of cosmological 
models on approach to spacetime singularities 

HI, 12], i. 

The idea that the fine structure constant, 
a, is a spacetime varying scalar field was first 
investigated by Bekenstein [3], who created a 
natural generalisation of Maxwell's equations to 
accommodate a varying electron charge. This 
idea was extended to include gravity and pro- 
vide a theory to explore cosmological conse- 
quences of varying a by Sandvik et al [5] . The 
resulting Bekenstein-Sandvik-Barrow-Magueijo 
(BSBM) isotropic cosmological models were 
found and used in conjunction with the astro- 
nomical data on varying a obtained from ob- 
servations of high redshift quasar spectra [6j. 
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More recently, the BSBM theory has been ex- 
tended to included the case where there is a 
coupling function (rather than simply a coupling 
constant) between the charged matter fields and 
the scalar field driving changes in a [7j and where 
that scalar field possesses a self-interaction po- 
tential 0. These theories are the analogues of 
the Jordan-Brans-Dicke theories for varying G 

In [10' it was shown how theories of this type 
could produce singularity-free homogeneous and 
isotropic cosmologies which displayed stable os- 
cillations around an Einstein static universe be- 
cause the effect of variations in the scalar field 
driving variations in a is to introduce a neg- 
ative 'ghost' density. Barrow and Tsagas [TT] 
considered a broader context for these solu- 
tions and showed how the inclusion of simple 
anisotropic expansion can modify the results be- 
cause the anisotropy can diverge just as quickly 
as a bounce-producing ghost scalar field on a 
approach to the singularity. In this paper we 
will consider more general closed anisotropic cos- 
mologies with anisotropic 3-curvature in this 
same context. 

Matter bounces introduced by the presence 
of ghost fields are not a new discovery (for a de- 
tailed examination see |l2j). However, in BSBM 
models the ghost field is an effective manifes- 
tation of underlying physics, not a new mat- 
ter source introduced by hand. In such mod- 
els quantum effects are ignored because of the 



prevailing attitude that ghost fields should not 
be quantized (and are in fact ill-behaved when 
quantized, with negative probability states). 
Furthermore, any coupling between a ghost field 
and a non-ghost field would allow an infinite 
amount of energy to be transferred from the 
ghost field. 

We simply take the view that BSBM models 
can serve as test models for bouncing cosmolo- 
gies. The idea of a "phoenix" universe within 
relativity is almost as old as big-bang models 
themselves, and goes back to Tolman [13] and 
Lemaitre IHj. This classical picture of oscillat- 
ing closed universes with zero value of the cos- 
mological constant, A, painted by Tolman is well 
known. If there is no entropy production then 
cycles for the time-evolution of the scale factor 
are periodic with the same amplitude and to- 
tal lifetime. If entropy increase is introduced 
in accord with the second law of thermodynam- 
ics then the oscillating cycles become larger and 
longer to the infinite future. The classical pic- 
ture for isotropic universes was competed by 
Barrow and Dcibrowski [15j, who showed that if 
a positive cosmological constant is included then 
the sequence of growing cycles will always come 
to an end, no matter how small the value of A > 
0. The ensuing behaviour will be to approach de 
Sitter expansion. If the entropy increase from 
cycle to cycle is small then the asymptotic state 
will be one in which the expansion is very close 
to a zero-curvature state with comparable en- 
ergy densities associated with matter and dark 
energy (ie the cosmological constant). The dark 
energy will necessarily be slightly dominant and 
the curvature will be positive - not unlike the 
situation in our observed universe. Barrow and 
Dcibrowski [TB] also considered the evolution of 
some simple bouncing anisotropic universes of 
Kantowski-Sachs type, but not in a context that 
included varying constants. 

Many current quantum theories of gravity ex- 
hibit curvature singularity avoidance, often in 
the form of a bounce (although they do not 
necessarily avoid geodesic incompleteness). In 
Loop Quantum Cosmology, holonomy correc- 
tions to the Friedmann equation give rise to a 
bounce at Planck scales (see p^ for a review). 
Horava-Lifschitz gravity also introduces higher- 



order curvature corrections to Einstein's equa- 
tions which can cause the universe to bounce 
[17J for some parameter choices. In the latter 
case the dynamics of an anisotropic solution have 
also been explored [18] . 

The aim of this paper is to extend [lOj and 
[15] to spatially homogeneous models which ex- 
hibit local rotational symmetry (LRS) an so 
are effectively axisymmetric. LRS models ex- 
hibit some of the features of full anisotropic 
model [19], yet the differential equations gov- 
erning their dynamics can be solved with rela- 
tive ease using numerical solvers and exact meth- 
ods. We will begin in section |ll] by setting our 
the action principle underlying variation of the 
fine structure 'constant', then in section IIIII we 



set out the equations of motion for our system 
and define quantities of physical interest, such 
as shear and Hubble expansion rates, in terms 
of metric variables. Section IIVI deals with two 
specific solutions to the equations of motion: a 
static solution and ghost-induced inflation. In 
particular, we will focus on the role played by 
anisotropies in both these cases and examine 
perturbations about isotropic cases. 



II. THEORIES OF VARYING ALPHA 

Varying "constants " can be described by 
extensions of the standard model of particle 
physics and/or general relativity (GR) by the 
promotion of constants to space and time de- 
pendent scalar fields. A well known and much 
studied example is that of Jordan-Brans-Dicke 
theory in which GR is extended by generalising 
Newton's constant G to become a field variable 
[9]. These self-consistent models for the varia- 
tion of constants necessarily contain conserva- 
tion equations for the energy and momentum 
carried by the varying scalar field and the grav- 
itational field equations account for the scalar 
field's effect on the spacetime geometry. This 
is in contrast to much of the old literature on 
varying constants, other than G, which merely 
'write-in' variations of constants into the equa- 
tions that hold in the theory where the constant 
does not vary. The existence of a self-consistent 
theory for the variation of a constant also shows 



that much discussion about the meaning of the 
variation of dimensional constants is not relevant 
because the solution of the second-order conser- 
vation equation for the scalar field describing the 
variation of a traditional constant always pro- 
duces constants of integration with the same di- 
mensions as the varying constant and a dimen- 
sionless combination is trivially available. 

Physical models with extra dimensions of- 
ten exhibit massless or light degrees of freedom 
which can lead to the variation of such constants 
|20j , [21j and there is the possibility for observa- 
tional bounds to be placed on any shift in the 
size of extra dimensions over the age of the uni- 
verse p2].p3]. 

In this paper we shall follow the BSBM model 
in which the fine structure constant, a is taken to 
be dynamical. Evidence of a dipolar spatial vari- 
ation has been recently claimed [6] and therefore 
it is natural to extend this scenario to consider 
space and time variations but in this paper we 
will only discuss time variations so that we can 
confine attention to ordinary differential equa- 
tions. Such variations are bounded by terrestrial 
experiments to have small variation at present 
|24j |25j . However evidence that the variation 
is small currently does not rule out more signifi- 
cant changes in the past. In particular, in BSBM 
theories a is not expected to vary during the ra- 
diation era, to increase only logarithmically in 
time during the cold dark matter dominated era, 
and then to become constant after the universal 
expansion begins accelerating. Thus laboratory 
experiments today would not be expected to find 
evidence for the variation of a found in high- 
redshift quasar observations (that derive from 
epochs before the universe began accelerating) 
even though it has been proved that any cosmo- 
logical variations in a will be seen in terrestrial 
experiments [26j . 

The BSBM model describes the effect of vary- 
ing the fine structure constant by the introduc- 
tion of a scalar dielectric field, ip with evolution 
of the charge of an electron given by e = eoe"^, 
in which cq is the value of the electron charge 
at some fixed time, for example today. Notice 
that Co is a fundamental constant and e/eo is di- 
mensionless. It has been shown [27J that in spite 
of modifying black hole solutions, the variability 



of a respects the second law of thermodynam- 
ics. This will be important in section III, as we 



will assume that all couplings between our fields 
obey the second law. 

The physical action is given by 

S= f V^{C, + Cm + L^ + e-^^Cem) (2.1) 

where Cg = R/WttG is the usual Einstein- 
Hilbert Lagrangian, i2^ = —^dutpd'^^ governs 
the scalar dielectric field, ip, Cem = —^f/iuf'^'^, 
and Cm is a matter Lagrangian independent of 
■0. Of particular importance to this paper is the 
constant coupling parameter u which we shall 
take to be negative, so rendering -0 an effective 
ghost scalar field. We do not consider the gener- 
alised case where UJ'{^p) 7^ 0, see |7j. From now 
on it will be convenient to simply consider the 
effective field, rather than the underlying dielec- 
tric. This field is massless, and it is clear from 
the action that its motion will be monotonic, as 
will be that of the effective induced fine structure 
constant. In terms of fiuids, this field will appear 
to be stiff, with equation of state Pip = p^ < 0. 
In a closed anisotropic cosmological model we 
will allow energy exchange to occur between the 
■0 field and an equilibrium radiation field with 
equation of state 3pr = Pr, to model an entropy 
increasing non-equilibrium process. 



III. COSMOLOGICAL EXPANSION 

The physical system under consideration will 
consist of a homogeneous anisotropic cosmology. 
For simplicity we will examine a system which 
is locally rotationally symmetric, and use this 
to gain insight into the more general case. For 
a concise review of these cosmological models, 
see J19j . This model is general enough to con- 
tain the purely general relativistic ingredient of 
anisotropic 3-curvature, which is missing from 
the simple anisotropic models of Bianchi types 
I and V. It includes the closed Bianchi type IX 
universe but only in the axisymmetric case where 
no chaotic behaviour occurs. The LRS type IX 
metric is 



where a* are the 5*0 (3) invariant 1-fornis |28j 



a = cos ipdO + sin if; sin 
a = — sin ^pd9 + cos tp sin 
a = dif) A- cos 6d(f) 

and the LRS condition requires 



These reduce to the special case of the 
closed Friedniann universes when a = b. The 
shear does not evolve with a oc (a6^)~^ as in 
Bianchi type I models because of the 3-curvature 
anisotropy on the right-hand side of eq. (3.4). 



For ease of exposition, let us define the variables 
r = ^ and Ht, = 6/6. The essential field equations 
then simplify to 



hij = diag{a{t),b{t),b{t)} 

The metric contains two time-dependent 
scale factors (due to the LRS condition), a{t) 
and 6(t). The energy densities are denoted by 
Pr for radiation, p^ for the scalar field, and p\ 
for the cosmological constant, and the total den- 
sity and pressure are p and p, where 



P = Pr + Pi, + PA- 

All matter matter sources have isotropic pres- 
sures. The independent variables are the princi- 
pal 3-curvatures 
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the mean Hubble expansion rate is defined by 
and the expansion shear scalar by 



a 



3^6 



These variables are subject to a constraint 
equation (the generalized Friedmann equation 
with SttG = c = 1) 
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^~62"464+^^ + 62- 
The remaining field equations are: 
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3a = Hh — 3a. 



(3.6) 



(3.7) 
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The continuity equation, which implies con- 
straint conservation, is 



p + 3H{p + p) = 0. 



(3.9) 



This governs the total energy density and 
pressure. Now, we introduce some energy ex- 
change between the fluids so that they obey 

p'^ + 6Hp^ = s, (3.10) 

Pr + 4.Hpr = -S, (3.11) 

where, s parametrises the flow of energy be- 
tween the scalar fleld and radiation. For our 
purposes, this will be taken to be of the form 



where 



-Pi>P 



/3 = /30 + PhH^ + Pa(T^ 



can include a linear coupling /3o, plus possible 
bulk, Ph, and shear, 13^, viscous contributions. 
In general, the /3's need not be constants. 

The scalar dielectric field evolves according 
to 



i^ + {3H + (3)iP = 



(3.12) 



and so 



tj) (X a exp[ 



pdt]. 



IV. SOLUTIONS 



In what follows we consider the Bianchi IX 
case with no cosmological constant (pA = 0)- If 
the matter content is a perfect fluid and obeys 
p + 3p > 0, then these universes expand from 
an initial curvature singularity to a maximum 
size before collapsing back to a future curva- 
ture singularity; the spacetime is past and future 
geodesically incomplete. However, since the ef- 
fect of varying the fine structure constant is to 
produce a ghost scalar field with p^ < which 
will dominate dynamics at small length scales, 
the energy condition is violated and solutions to 
the BSBM model exist which have infinite past 
and future temporal range. 

In this section we will examine two particular 
solutions to the equations of motion. The first 
is that of a static spacetime, and its behaviour 
under perturbations. The second is that of an 
spacetime in which the coupling between fields 
leads to infiationary behaviour. 

A. The Static Solution 

There exists a static solution of the form 



e{t) 



eo sm 



(\/2t). 



(4.3) 
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(4.1) 



for any given value of 6. Note that in order to 
be static the solution must be isotropic {a = 0), 



since from eq. (3.4) the 3-curvatures must match 
{R\ = R2) for the shear to remain constant and 
eq. (|3.8l) requires o" = 0. 



Now consider making a small perturbation 
about the isotropic solution by introducing a 
small anisotropy: r = 1 + e where r = 1 rep- 



resents isotropy. Prom 3.6, we find that 



-3eHb + 



2e + 3e2 + e^ 
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(4.2) 



Without loss of generality, we take the un- 
perturbed static solution to have 6 = 1. If we 
introduce small parameters 5{t) and T]{t) such 
that 6 = 1 + rj, and p = 3/2 + 6, then to first 
order in our small parameters: 
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FIG. 1: e versus time with no coupling (s = 0). Ini- 
tial values: r — 1.05, a — 0, H — 0. 

In the case where there is no coupling be- 
tween the fields (s = 0), these oscillations con- 
tinue endlessly as shown in figure [TJ However, 
once coupling is introduced (s / 0), the static 
case becomes unstable because the balance be- 
tween p^ and pr is broken, energy is transferred 
from the ghost field that supports stable oscilla- 
tions and eventually, after several oscillations, 
it settles into radiation-dominated expansion, 
shown in figure [2| 
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FIG. 2: e versus time with coupling turned on, s 7^ 0, 
showing the system isotropizing. Initial values: r = 
1.05, cr = 0, F = 0, /3o = 0.05, Ph = /3a = 0. 



From the constraint equation 3.1, we find 
that to first order, our small parameters are re- 
lated by 26 = —e — 3ry. Furthermore, we can 



decompose 5 into the radiation and scalar field 
components, 5r and 6^. In the absence of field 
couplings, there is a relationship between these 
fields, due to their coupled equations of motion 



3.12 3.10 



3/2 

Pr OC p^ . 



(4.4) 



In the static case under consideration, the 
constant of proportionality is — v6. Note that 
this relationship is broken by introducing a cou- 
pling between the fields. For small 6 we are 
therefore led to Sr = 4:6 = — 2e — 6rj. 



From 3.3, the evolution of r] is given by: 
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Therefore there is a (more complicated) sta- 
ble oscillatory behaviour for tj about the static 
solution and the evolution of e has already been 



determined by 4.3 Thus we have an unusual be- 
haviour characterised by stable anisotropic oscil- 
lations around the isotropic Einstein static uni- 
verse. This generalises the simple isotropic os- 
cillations about the static universe that exist in 
Friedmann universes with a ghost field found in 
dOj and lllj. 

We can now determine further effects of al- 
lowing a coupling between the fields. The sec- 
ond law of thermodynamics requires s > 0. The 
exact form of s - taking into account terms rep- 
resenting constant coupling, bulk and shear vis- 
cosities, will of course affect the exact dynamics. 
However, it is possible to make progress by as- 
suming only that s is non- negative. 

First consider the case of no coupling (s = 0). 
The evolution of the radiation field is determined 
by 
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Pr 
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(4.6) 



Due to the cyclic behaviour of e and r/ (and 
hence of their derivatives), 5r will also cycle, re- 
turning to its initial value, as any integral of the 



right-hand side of 4.6 across a complete cycle will 
be zero. However, with a positive coupling be- 
tween the fields, this relationship is broken, and 




Time 

FIG. 3: Spatial volume versus time with coupling 
turned on (s ^ 0). The system oscillates, but ex- 
pands from one cycle to the next. Initial values: 
r = 1.05, cr = 0, iJ = 0, /3o = 0.01, (3h ^ Pa ^ 0. 



a term which is always non-negative (and so has 
a positive integral across cycles) must be added. 
Hence, across a cycle in 77 and e, the value of 5r 
now increases and we must adjust our equation 
of motion 14.51 to include this term. We write 



f/ = e + r/ + 



A, 



(4.7) 



where A is a positive term representing the in- 
crease in 5r due to field couplings created by in- 
troducing s > 0. The variable r] no longer cy- 
cles about zero, and the system is slowly pushed 
away from stability, and enters a pseudo-cyclic 
phase in which a series of bounces occur with 
increasing local minima and maxima of the ex- 
pansion volume, a6^, as shown in figure [3J The 



behaviour of e 4.2 is affected by this expansion 
(recall that e is already small). The expansion 
of h{t) means that Hi, is no longer small, and 
the equation gains a damping term. Similarly, 
the frequency of oscillations, \/2/6, is reduced by 
this expansion in h{t) and the solution takes the 
approximate form of a damped harmonic oscil- 



lator. Note that in the derivation of 4.2 only the 
smallness of e was used - therefore this damping 
behaviour is present in all expanding solutions. 
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FIG. 4: Hubble expansion rate versus time for an 
inflating solution. Initial values: /3o = 0.03, /3h = 
/3^ = 0, r = 1, CT = 5 X 10-4^ a = 10^, H = l.6x lO^^. 



B. Ghost-induced Inflation 

Spacetimes which exhibit inflation are of spe- 
cial interest to cosmologists because inflation can 
solve a number of well known puzzles about the 
universe's structure [29], and make a series of 
detailed predictions that can be tested by ob- 
servations of the microwave background radia- 
tion |30j . Typical inflationary models exhibit 
expansion in which the Hubble parameter is (ap- 
proximately) constant for a finite time interval. 
In general relativistic cosmology this is usually 
achieved by introducing a matter content that is 
(or is equivalent to) a scalar field subject to a 
self-interaction potential whose contribution to 
the total energy density is dominant during this 
expansion, with an equation of state close to that 
of that produced by an exact cosmological con- 
stant with PA = —PA- 

Inflation induced by ghost fields has been 
studied as an alternative to the usual slow-roll 
models [31]. Such models have potentially ob- 
servable consequences for the microwave back- 
ground trispectrum [32j [33j , but require that the 
translation invariance of the scalar ghost field 
is broken. In the BSBM models under consid- 
eration, however, translation invariance can be 
preserved, with the field coupling responsible for 
creating the inflationary energy density. 

Let us examine the case of a linear coupling- 
induced inflation with s = /3op^, for constant 
/3o > 0. When the volume is large there exists 
an asymptotic solution of the form 
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in the isotropic case. This solution is stable, 
and is approached by dynamical trajectories, as 
shown in figure |4j Under a small perturbation 
H = — /3o/2 -|- h and r = 1 -|- e we find that to 
first order in the small parameters: 
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(4.9) 
(4.10) 



In this solution, b{t) is exponentially growing, 
so the final term in |4.9| quickly becomes negligi- 
ble. Therefore, although a shearing expansion 
may occur, e quickly falls to zero, locking the 
shear at a fixed value. This infiationary phe- 
nomenon is not unique to linear couplings with 
Ph = /3o- = 0, but is simplest to demonstrate 
in this case. Likewise, there is no requirement 
for the matter field to consist solely of radia- 
tion - introducing more matter fields with cou- 
plings whose sign is determined to be in accor- 
dance with the second law of thermodynamics 
yields a system which also exhibits inflation of 
this type. For most couplings, however, the in- 
flationary phase will end when the dust fleld be- 
comes dominant. 

The evolution of the fine structure constant 
is shown in figure [SJ Initially, the solution is 
like an ascending staircase with rapid changes 
at each scale factor bounce, see for compari- 
son y^ . Monotonicity of a is ensured since the 
scalar field cannot have positive energy density; 

since ip = y ^* , we have ■(/' > for all time, 
and so ^ cannot oscillate through maxima and 
minima. Across repeated bounces, log{a) will 
appear to increase in steps when \p^\ is small 
as the relative size of this energy density os- 
cillates greatly within a single cycle. However, 
as energy is transferred into the radiation field, 
these steps will become less apparent, eventually 
approaching a constant gradient once the ghost 
field reaches the condition for de Sitter infiation 
to occur. Thus, even though the universe oscil- 
lates from cycle to cycle, the fine structure 'con- 
stant' continues increasing from cycle to cycle 
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FIG. 5: Evolution of log(a) versus time shown over 
a short timescale (top), and the same evolution over 
a long timescale (bottom). Initial values: /3o = 0.01, 
Ph = Pa = 0, r = 1.01,0- = 0,a = 1,H = 0,p^, = 
-0.01. 



and there will typically only be a finite inter- 
val of cycles in which a takes values that allow 
stable atoms to exist [Ml , ES] . 



V. DISCUSSION 

In this paper, we examined the new effect of 
introducing anisotropics into the BSBM frame- 
work for varying a, although the conclusions 
have broader applicability to anisotropic cos- 
mologies containing ghost fields and entropy- 
increasing energy exchanges between fields. In 
particular, we studied the dynamics of locally 
rotationally symmetric Bianchi IX cosmologies. 
It was shown that under certain conditions the 
bouncing behaviour observed in isotropic models 
persists, with the fine structure constant" chang- 
ing in an almost step-like increasing manner be- 
tween cycles as time increases. It is apparent 
from 13.21 that on short scales there is a tension 



between the shear terms and ghost field, as both 
scale as the inverse square of the volume. When 
the anisotropy is small, the contribution from 
the ghost field dominates. This leads to a bounc- 
ing model, reproducing closely the results seen in 
|10| . Furthermore, there exists a static solution, 
perturbations about which lead to a sequence 
of anisotropic bouncing phases. When there is 
a coupling between the matter fields, the second 
law of thermodynamics ensures that this process 
isotropizes the system by energy exchange. 

The resulting dynamics lead to a pseudo- 
cyclic universe in which the fine structure con- 
stant monotonically increases across bounces, 
and for small values of the associated dielectric 
scalar field, this increase is dominated by dynam- 
ics near the bounce point. The minimum and 
maximum volumes of the universe also increase 
across cycles, with the total energy density de- 
creasing. Eventually, the model reaches a point 
at which the coupling between fields fixes the 
energy density to be constant in time, and the 
universe undergoes a de Sitter phase in which it 
inflates. Since this model is limited to include 
only the dielectric field and radiation, there is 
no transition to the dust-dominated era that one 
would expect at the end of this phase, and so in- 
flation is endless within the model. It is possible 
to flnd solutions with dust in which the system 
again reaches a point of steady inflation. How- 
ever within the space of couplings with s > 
which obey the second law, these solutions are a 
set of measure zero. 

The evolution of a throughout the history of 
BSBM universes displays interesting traits. At 
late times, in a large universe it will appear that 
a has settled to a constant value. In doing so, 
throughout a series of oscillations the universe 
will have isotropised greatly, with a stepping up 
between cycles. As the dynamics are invariant 
under changing the initial value of a there is no 
obvious mechanism to determine the constant to 
which it will approach. 
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